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Abstract 

Multifractal analysis aims to characterize signals, functions, images or fields, via the fluctuations of their 
local regularity along time or space, hence capturing crucial features of their temporal/spatial dynamics. 
Multifractal analysis is becoming a standard tool in signal and image processing, and is nowadays widely 
used in numerous applications of different natures. Its common formulation relies on the measure of local 
regularity via the Holder exponent, by nature restricted to positive values, and thus to locally bounded 
functions or signals. It is here proposed to base the quantification of local regularity on p-exponents, a novel 
local regularity measure potentially taking negative values. First, the theoretical properties of p-exponents 
are studied in detail. Second, wavelet-based multiscale quantities, the p-leaders, are constructed and shown 
to permit accurate practical estimation of p-exponents. Exploiting the potential dependence with p, it is also 
shown how the collection of p-exponents enriches the classihcation of locally singular behaviors in functions, 
signals or images. The present contribution is complemented by a companion article developing the p-leader 
based multifractal formalism associated to p-exponents. 

Keywords: Pointwise regularity, p-exponent, wavelet p-leaders, negative regularity, singularity 
classihcation, multifractal analysis 


1. Introduction 

Context: multifractal analysis, local regularity and Holder exponent. Multifractal analysis 
nowadays constitutes a classical signal and image processing tool, available in most modern toolboxes. It 
is commonly used to model or analyze scaling properties and local irregularities in signals and in image 
textures. It has been involved in a large variety of real-world applications of very different natures, ranging 
from biomedical (heart rate variability [1], neurosciences [2, 3]), to physics (turbulence [4]), geophysics 
(rainfalls [5], wind [6], earthquakes [7]), hnance [8, 9, 10], music [11], or Internet traffic [12], to name but a 
few. 

In essence, multifractal analysis describes, via the multifractal spectrum, the huctuations of local regular¬ 
ity along time, or space, of a function, signal, image or held. Fundamentally, it relies on two key ingredients: 
At the theoretical level, a pointwise regularity exponent h{t) that formalizes the intuition of local regularity; 
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at the practical level, a multiscale quantity Tx{a, t) that permits the actual measurement of local regularity, 
via a local power-law behavior as a function of the analysis scales a: 

Tx{a,t) when a —0, 

or, more technically, 

a-s-0 logo 

In its commonly, if not exclusively, used formulation, multifractal analysis relies on the quantification of 
local regularity via the so-called Holder exponent [f3, 14, 15, 16]. In practice, several multiscale quantities 
have been involved in multifractal formalisms, the practical counterpart of multifractal analysis. It has 
now been long recognized that increments or wavelet coefficients do not yield accurate analysis of Holder 
exponents. The earliest solid practical formulation of multifractal analysis was based on the (continuous) 
wavelet transform modulus maxima method (WTMM) (cf. e.g., [17, 18, 19, 20]). A multifractal extension 
of Detrended Fluctuation Analysis (MF-DFA) has also been proposed (cf. [3] for the founding article, and 
[21, 22] for related developments). Recently, a theoretically well-grounded and practically efficient formula¬ 
tion that, unlike WTMM and MFDFA, extends well to higher dimensional signals, has been proposed: It 
relies on wavelet leaders, constructed as local suprema of discrete wavelet coefficients, cf. [16, 23, 24, 25]. 

Limitation: negative regularity. An important, yet often overlooked, issue in multifractal analysis, 
consists in the fact that choosing multiscale quantities on which analysis is based also amounts to selecting a 
specific measure of pointwise regularity. For instance, the wavelet leaders are the only multiscale quantities 
that were shown theoretically to be able to actually characterize the Holder exponent [16, 23, 24, 25]. 
Conversely, the practical use of a particular regularity exponent necessitates the construction of multiscale 
quantities specifically tailored to it, requiring the verification of global a priori regularity assumptions, which 
may not always hold for the data to be analyzed. In particular, by definition, Holder exponents cannot take 
negative values. In its current formulation, multifractal analysis is thus restricted to functions whose local 
regularity is everywhere positive, i.e., to functions that are everywhere locally bounded and thus show a 
minimal global regularity. This is a severe restriction for the practical use of multifractal analysis as signals 
and images from real-world applications are often characterized by discontinuities and thus by a negative 
uniform, or minimum global, regularity. This is notably quasi-systematically the case for biomedical time 
series, (cf. e.g., [26, 27] and references therein for reviews). 

It has often been proposed to circumvent this well-recognized limitation by integrating data prior to 
performing multifractal analysis [17, 18, 19, 20]). The use of fractional integration, of order tuned to the 
uniform regularity of data, has recently been proposed and embedded directly in the definition of wavelet 
leaders [25, 26, 27]. However, fractional integration yields the significant issue of relating the characterization 
of local regularity of integrated functions to those of the original functions, prior to integration. Notably, 
the commonly used rule of thumb that fractional integration results in a uniform shift of the multifractal 
spectrum is in general wrong (specifically, in the presence of oscillating singularities [25]). This issue will be 
further discussed in Section 4. 

Goals, contributions and outline. In this context, as a possible alternative to (fractional) integration, 
we propose the use of p-exponents, which potentially take negative values and hence permit to characterize 
negative local regularity. Though introduced in the theoretical context of PDFs as early as 1961 for p > 1 
by Calderon and Zygmund [28], p-exponents were not used in signal processing until the 2000s when their 
wavelet characterization was proposed [29, 30, 31]. In the present contribution, we study which information 
on the local behavior of a function near a singularity is supplied by the knowledge of the collection of 
p-exponents. 

After recalling definitions, Section 2 studies the theoretical properties of p-exponents. 

The corresponding multiscale quantities, the p-leaders, are defined in Section 3 as local P norms of 
(discrete) wavelet transform coefficients. They will be shown to permit an accurate practical measure of 
p-exponents as well as to extend the definition of p-exponents to p S (0,1). Several representative examples 
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will be used to illustrate the properties and potential of p-leaders and p-exponents to characterize local 
singular behaviors. 

The joint use of all p-exponents paves the way to a thorough classification of local singular behaviors 
(in terms of p-invariant, canonical, oscillating, balanced or lacunary singularities); this will be detailed in 
Section 4. The interplay with negative regularity will also be discussed (Section 4.4). 

While the present contribution addresses pointwise or local regularity and hence focuses on isolated 
singularities, the final goal is multifractal analysis, aiming to study collections of intertwined singularities. 
Therefore, the present work is complemented by a companion article which develops and studies the cor¬ 
responding p-leader based multifractal formalism that permits the practical implementation of multifractal 
analysis based on p-exponents. 

Matlab routines designed by ourselves and implementing p-leaders and estimation procedures for p- 
exponents will be made publicly available to the research community at the time of publication of the 
present article. 

2. Pointwise regularity exponents 

Throughout the present article, {X{x),x e K'^} denotes the function or sample path of a stochastic 
process, or random field, in dimension d, to analyze. 

2.1. Pointwise Holder regularity 

The most commonly used notion of pointwise regularity is defined via the so-called Holder exponent, 
whose definition and properties are recalled below. For Section 2.1, {^(a;)}a;gR<i is assumed to consist in a 
locally bounded function. 

Local Holder spaces C°‘{xo). Let o > 0. The function X is said to belong to C°‘{xo) at location 
xq G with a > 0, if there exist a constant C > 0 and a polynomial P^g of degree less than a, such that, 
for a small enough, 

\X{xo + a) - Pxoixo + a)\ < Clal^. (3) 

Because it obviously generalizes the Taylor polynomial for functions, Px^ is usually also referred to as 
the Taylor polynomial of order a oi X at Xg. When a < 1, the Taylor polynomial boils down to a constant 
Pxo{x) = X(xo). A more general discussion on Px„ will be given in Section 2.2 (see also [32]). 

Pointwise Holder exponent. The Holder exponent of X at location Xq is defined as 

h{xo) = sup{a : X G C°'{xq)}. (4) 

It is commonly used to characterize the local regularity of X at Xg: The larger /i(xg) the smoother X around 

Xg. 

Limitation of pointwise Holder exponent: positive regularity. By definition, the Holder exponent 
/i(xg) cannot take negative values and (3) implies that X is bounded in a neighborhood of xg. Indeed, 
Pxg is bounded and, since a > 0, so is |a|“. However, large classes of signals and images from real-world 
applications cannot be satisfactorily modeled by locally bounded functions (cf., e.g., [27] where a practical 
wavelet based criterium to assess local boundedness is supplied and applied to numerous real-world data). 

The attempt to characterize singularities with negative regularity by simply allowing a in (3) to take 
negative values does not yield a satisfactory definition: Indeed, if (3) is satisfied for a given a (which may 
be negative), then, in any corona 0 < Ci < |x — xg| < C 2 , / is a bounded function. Therefore, such 
a definition would only permit to define isolated singularities of negative order, which would thus not be 
relevant in the framework of multifractal analysis where singularities with a given exponent may be dense, 
and usually consist of sets with strictly positive Hausdorff dimensions (whereas a set of isolated points is at 
most countable and hence with Hausdorff dimensions equal to 0). 
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2.2. p-exponent regularity 

To characterize negative regularity in data, a new definition for pointwise regularity is used, the Tp^xq) 
regularity, which has the advantage of relying on the less restrictive assumption that data locally belong to 
LP{R‘^), instead of requiring local boundedness. This notion was introduced by A. Calderon and A. Zygmund 
[28] and has recently been put forward in the mathematical literature in [29] (see also [31, 26]). It permits 
the definition of a collection of p-exponents to measure pointwise regularity. 

Definition 1. Let p > 1 and X be a function that locally belongs to Let B{xo,a) denote the ball 

centered at xq and of radius a > 0. Let a > —d/p. A function X belongs to TP{xo) if there exist a constant 
C and a polynomial Px„ of degree less than a such that, for a small enough, 

T^P\a,Xo) := (\ [ |A(u + xq) - <Ca°‘. (5) 

y a JB{0,a) J 

The p-exponent of X at Xq is defined in [29] by: 

hp{xo) = supja : X e Tf/fx^)}. (6) 

The p-exponent hp{xQ) takes values in [—d/p, -l-oo] (see Theorem 1 below), thus allowing for a proper math¬ 
ematical dehnition of negative order pointwise regularity. This definition constitutes a natural substitute 
for the pointwise Holder regularity when dealing with functions which are not locally bounded. The usual 
Holder regularity actually corresponds to p = -l-oo (in this case, the local L^ norm in the left hand side of 
(5) boils down to a local norm, thus supplying a condition equivalent to (3)). Definition 1 does not 
include p-exponents with p < 1 because L^ with 0 < p < 1 is mathematically ill-defined. Section 3.4 will 
however show how the wavelet framework permits to extend the definition of p-exponents to 0 < p < 1. 


Taylor polynomial. The polynomial Pxq defined implicitly by (5) is unique for a given a independently 
of p. When a crosses an integer value N , P^g may be modified by addition of terms of degree N. However, 
for two different Oi and a^-, the expansions of P^g coincide up to order min(Q;i,a 2 ) (cf- Appendix A). This 
implies that by picking up the integer part of hpix^) we get the polynomial P which corresponds to the 
largest possible value of a. Theorem 1 (Point 2) below indicates that one can fix a unique Taylor polynomial 
of X at Xo, whose coefficients are independent of p, and are referred to as the (generalized) Peano derivatives 
of X at xq [33, 32]. 

One of the main advantages of the wavelet framework developed in Section 3 is however that the compu¬ 
tation of P^g is not required to measure the p-exponent. The Taylor polynomial is thus not further discussed 
and readers are referred to e.g., [32]. 


2.3. Properties of p-exponents 

Let us now state the main theorem characterizing the properties of the collection of p-exponents, or 
mapping p —>• hp{xo), at location xq. Let us also introduce the notation X G LP{R‘^) to indicate that X 
belongs to LP{M.‘^) locally. 

Theorem 1. Let X G L)^j,(M‘^) and xq G K'^. 

Let po{xo) = supjp : X G L/^^{M.‘^) in a neighborhood of xq\. 


The function p —>■ hp(xo) is defined on [l,po(a;o)) (and possibly also at po(xo)), and possesses the following 
properties: 


1 . 

2 . 

3. 


It takes values in 


d 

Poi^o) ’ 


OO 


It is a decreasing function of p. 

The function r i—)■ hi/T^xe/) defined on [0,1] 


is concave. 
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Furthermore, Conditions 1 to 3 are optimal, i.e. if po G (l,c») and (j) is a function defined on [l,po] 
satisfying the above conditions, then there exists X G LP° such that 

Vp e [l,Po], hp{xo) = (j){p). (7) 

The proof of the first part of Theorem 1 is detailed in Appendix A. The proof of optimality is detailed 
in Appendix B and relies on the explicit construction of a function X such that the function p —>■ hp{xo) 
can be any function satisfying Conditions 1 to 3 from Theorem 1. This explicit construction also serves 
as a reference example of singularities where p-exponents differ at Xg = 0 (i.e., the function hp(0) is not a 
constant), see Section 2.4.3 below. 


Theorem 1 calls for the following remarks: 

First, as already stated above, p-exponents can by definition take values down to —d/p and hence allow 
to formalize the notion of negative regularity exponents. 

Second, p-exponents for different values of p do in general not coincide. In particular, the Holder exponent 
(p = oo) does not in general coincide with p-exponents for p < oo. 

Third, as a consequence of the concavity of the function r i—>■ hi/r{xo), it follows that, when it takes 
finite values, hp(xo) as a function of p is continuous, except perhaps at endpoints 1 and po (indeed a concave 
function is continuous except perhaps at the end-points of its domain of definition). 

Finally, the more difficult problem of understanding which properties are satisfied by the functions of 
two variables (xo,p) -G hp(xo) remains largely open. 

2 . 4 . Pedagogical and reference examples 

Let us now illustrate p-exponents on pedagogical examples. 


2.4.I. Example 1: Cusp 

The very reference for local singularity consists of the cusp function: 

If a ^ 2N Cq(x) = |a; — xo|“ (8) 

It is straightforward to show that the p-exponent of Ca at Xq does not depend on p and takes the constant 
value a, as is the case with Holder exponent, i.e., 


hp(xo) = a, 


for all p < Po with 

- 1-00 when a > 0, 

—d/a when a < 0. 

Examples of cusp singularities with positive and negative regularity are illustrated in Fig. 2. 


Po = 


( 9 ) 


2 . 4 . 2 . Example 2: Chirp 

Chirp functions are another classical and commonly studied example of singular behavior. They are 
defined as 

Ca,p{x) = \x-XQ\°'sm.(- -^— 73 ) ( 10 ) 

\\x-xo\Pj 

and serve as the reference for the class of oscillating functions. The application of (3) immediately yields 
that 

hp{xo) = a, 


for all p < Po with 


Po 


-|-oo when a > 0, 
—d/a when a < 0. 


( 11 ) 


Two chirps with positive and negative exponents are illustrated and analyzed in Fig. 3 (cf. Section 4.3 for 
further discussions). 
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2.4-3. Example 3: Lacunary comb 

For cusp or chirp singularities, p-exponents do not vary with p. To illustrate the importance of the 
possibility of making use of a collection of different p-exponents, let us now introduce constructive examples 
of univariate functions ^^,7 : K —>■ K, referred to as lacunary combs, which highlight a major benefit of 
using p-exponents: They enable the refined characterization of certain types of singularities that cannot be 
revealed when using the Holder exponent. These functions constitute a key-construction in the proof of 
Theorem 1. 

Let a G M and 7 > 1. The function Fa^-y, defined as: 

_ r2-“J for a; G [2-^,2-^ -t 2-^i] for j G N+, 

10 otherwise. 


has support on a set U.y = U^>p[2“-^, 2“-^ -|- Examples of functions with a > 0 and a < 0 are 

plotted and analyzed in Fig. 4. 

The goal is now to derive the function p —>■ hp{0) for Fa^-y at location xq = 0 and to compare it with the 
Holder exponent. 

When a > 0, Fa^j is bounded in the neighborhood of xq — 0, so that it belongs locally to all spaces 
with p G [l,-|-oo]. The function Fa^j is continuous at 0 where it vanishes (see Appendix B). Moreover, 
|^'ct. 7 (a^)| < |a:|“, and this estimate is an equality at points 2~F Therefore, the Holder exponent reads 
h{0) = a. 

When a < 0, let a > 0 and J such that 2“*^ -|- 2~^'^ < a < 2“'^+^; then 


1 

a 


\Fa,j{x)\Pdx = ^ V2 
a 


--ijry-apj 


(13) 


For any p < the sum (13) converges and therefore Fa^-y belongs locally to F^ . The function Fa^-y is no 

longer bounded near xq = 0, and the Holder exponent at xq = 0 is thus no longer defined, yet p-exponents 
are well-defined as long as p < — 7 / 0 :. Note that the minimal regularity condition / G requires that 
a > — 7 , which we assume from now on. 

Let us compute the p-exponents of F^^y at Xq = 0, assuming that p < —^ja if a < 0. The relation (13), 
with 2~’^ + "£ a < 2“'^+^, implies that 

- / |F„ JxWdx = - V 2 -'^^2-“P^ ~ 

a J-n ' a 


This yields 

T^^^(a,xo) ~ (14) 

If we now take a in the interval [2~^, 2“-' -|- 2“'''-^], we obtain a quantity which is bounded by the value of 
T^\a,Xo) in (5) at 2“1 -|- 2“'''!, so that (14) still holds. We have therefore obtained that, Va > — 7 , the 
p-exponent of Fa^j at the origin is 


with 


Vp G [l,po) 


hp{0) — a. - 

P 


Po 


- 1-00 when a > 0 , 

—^jcx when a < 0 . 


(15) 


(16) 


The theoretical values and practical estimates of hp{0) are illustrated in Fig. 4. 

This example shows that the p-exponents of Fa^j differ at xq = 0 for all values of p G [l,Po) and thus 
supply a much richer characterization tool, even in the case when a > 0 and the Holder exponent can be 
computed. This will be further discussed in Section 4.3. 

Here, we have constructed a function such that the p-exponent is an affine function of 1/p; the general 
case (pertaining to the proof to the second part of Theorem 1) is detailed in Appendix B. 
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3. Wavelet characterization 


3.1. Discrete wavelet coefficients and p-leaders 

It has long been recognized that wavelet coefficients constitute ideal quantities to study the regularity of 
data (see e.g., [34, 18, 14, 15, 16]). The characterization of p-exponents proposed here relies on the use of 
the d—dimensional discrete wavelet transform (dDWT), whose definition and properties are briefly recalled. 

Mother wavelets. Let denote a family of oscillating functions with fast decay and 

strong joint time-frequency localization properties, referred to as the mother wavelets. Let us further assume 
that these functions are chosen such that the collection of templates of 

2 -47/2^W(2-Ja: _ k), for i = 1, • • • 2^ - 1, j gZ and kG (17) 

dilated to scales a = 2-1 and translated to space positions 2^k, form an orthonormal basis of [35]. 

The d-dimensional orthonormal wavelet bases mostly used in practice are obtained by tensor product of 
univariate orthonormal wavelet basis [36]. 


Discrete Wavelet Transform. The coefficients of the dDWT of X are defined as 


cfl = f X{x) {2-^x - k) dx. 


(18) 


Note the use of an ffi normalization for the wavelet coefficients that better fits local regularity analy¬ 
sis and yields the correct self-similarity exponent of the wavelet coefficients for self-similar functions, see 
[18, 16, 24, 32]. For further details on wavelet bases and wavelet transforms, the reader is referred to e.g., 
[37]. 


Uniform regularity and number of vanishing moments. The mother wavelets {X)}i=^,...2d-l 
are further required to possess additional regularity and localization properties: They are assumed to belong 
to with ry, as large as possible. When r^p > 1, all their partial derivatives of order at most have 

fast decay. Also, the number of vanishing moments Np, is defined as a positive integer such that for any 
polynomial P of degree strictly smaller than Np,., 


P{x)ip^'^\x)dx = 0 . 


(19) 


Both the regularity and the vanishing moment assumptions are required in order to obtain accurate wavelet 
characterizations of pointwise regularity: Let /imax denote the largest smoothness order found in X, then a 
sufficient condition for choosing the mother wavelet reads: 

^max ^ min (^rpi^N^p). (^0) 

In general, one does not have information concerning a priori regularity of the data, and therefore, one 
does not know how smooth the analyzing wavelets should be. In practice, a rule of thumb consists in 
using smoother and smoother wavelets, until the outcome no longer depends on the wavelet used, which 
is interpreted as meaning that sufficient regularity has been reached. This can afterwards be confirmed 
using multifractal analysis tools (see [38]), which allow to determine the maximum regularity exponent 
present in the data. Further, with orthornormal wavelet bases (such as the so-called ’’Daubechies” compactly 
supported wavelets [39] , widely used in applications, and used in this contribution as well as in the companion 
article [38]), one necessarily has N^j > r^. A sufficient (and conservative) condition for accurate wavelet 
characterizations of pointwise regularity simplifies to /imax < rp,. 

Another (related) consequence of the vanishing moment requirement of major practical importance con¬ 
sists of the fact that wavelets are orthogonal to polynomials: The wavelet characterization of either Holder 
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Figure 1: Definition of p-leaders. 


or TP{xo) regularity exponents thus avoids the explicit estimation of the Taylor polynomial. 


p-leaders. The multiscale quantities suited to the characterization of p-exponents are referred to as p- 
leaders and constructed as the norm of a subset of wavelet coefficients c\. Let k = (fci,..., kd) and 
j G 'Z, denote a dyadic cube 

A,,fe := 2^ (fci + 1)) X • • • X [2^kd, V{k^ + 1)) . 

Let CX denote the cube homothetical to the dyadic cube A, with same center, but C times wider. Accord¬ 
ingly, CA = cpk, and V'a = V'yl- 


If X e the p-leaders are defined as: 





i/p 


2-dU-3') 


( 21 ) 


where j' < j is the scale associated with the sub-cube A' included in 3A (i.e. A' has width 2-i ). The £p 
norm is thus taken on all cubes A' of scale at most 2^, which are included either in A or in its — 1 closest 
neighbors, as illustrated in Fig. 1. 

When p = - 1-00 (and X G L“^(IR‘^)), p-leaders boil down to the classical wavelet leaders 


h = 




sup I c\ 

j'<j, A'C3A 


used to characterize Holder exponents [16, 24]. 

Remark. Note that an equivalent definition of p-leaders is given by 


&\j,k) 



i/p 




( 22 ) 


where dj^k = supj_]^ 





3.2. Global regularity and wavelet coefficients 

Local boundedness. As already discussed, the Holder exponent suffers from a well-known and docu¬ 
mented limitation: It can be used as a measurement of pointwise regularity for locally bounded functions 








only, see [40]. Whether empirical data can be well-modeled by locally bounded functions, or not, can be 
practically determined through the computation of the uniform Holder exponent /i™™, using the following 
wavelet characterization: 


log 


= liniinf 

j-f-oo 


sup 

k 




log(2f) 


(23) 


Indeed, if h""" > 0, then X is a continuous function, whereas, if < 0, then X ^ see [27, 26]. For 
numerous real world applications, the restriction > 0 constitutes a severe limitation, cf. [27]. From a 
practical point of view, mother wavelets satisfying r^ > are required for an accurate estimation of h™™. 


Local LP regularity. Wavelet coefficients can furthermore be used to assess whether X locally belongs 
to LP or not, i.e., to verify the condition that is a priori required for using the corresponding p-exponent 
[26, 32, 27]. Let Sc{j,p) denote the wavelet structure function, defined as the space/time averages of the 
magnitude of wavelet coefficients raised to a positive power p > 0: 


S.{j,p) = 2^^Y.E\c%\^ 

k i—1 


Let rjx (p) denote the wavelet scaling function, defined as 


Vp > 0, 


Vx (p) = lim inf 

j-f-oo 


log {Sa{j,p)) 
log(2f) 


(24) 


(25) 


Because of its function space interpretation in terms of Besov spaces (see Section 3.4 below), t]x{p) does not 
depend on the wavelet basis as long as > | | and p > 1. Then when p > 1: 


If Vx{p) > 0 then X e 1 

if Vxip) < 0 then X ^ j 


The wavelet scaling function rjx {p) must therefore be computed prior to applying p-exponent character¬ 
ization to determine the range of values of p suitable for analysis. 


3.3. p-exponent regularity characterization with p-leaders 

Let p > 1. When r]x{p) > 0 and r.,p > hp(xo), the p-exponent hp{xo) defined in (6) can be recovered 
from p-leaders [31, 26, 29]: 


hp{xo) = lim inf 

j-f-OO 


log (41(xo)) 

log(2f) 


(27) 


Note that the characterization provided by (27) does not require the computation of the Taylor polynomial, 
and that it extends to any p the previously obtained wavelet leader characterization of Holder exponents 
(valid when > 0) [16, 24]: 


h(xo) = lim inf 

j-f-OO 


log (^A,.,(.o)) 
log(2^) 


3 . 4 . The case 0 < p < 1 

Limitation p > 1. The restriction p > 1 may constitute a severe drawback in applications as it implies 
that pointwise singularities with regularity smaller than —d cannot be dealt with. As already pointed out, 
the restriction p > 1 stems from spaces being ill-defined when 0 < p < 1. In particular, even when 
the wavelet basis belongs to the Schwartz class and X satisfies / \X{x)\p dx < 00 ior p < 1, its wavelet 
coefficients might not be well-defined: Indeed for example, for the function X{x) = which belongs e.g. 
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to the scalar product f Xip is in general nndefined even if ip is C°°, unless (/^(O) = = 0. 

Replacing Lp spaces with Besov spaces. However, this issue has been dealt with mathematically by 
the substitution of spaces by Besov spaces, see [36, 35]. We sketch below only the basic notions together 
with practical implications that are relevant to pointwise regularity. 

Let p > 0 and s G K. A function X belongs to the Besov space if its wavelet coefficients satisfy 

3(7 Vj 2^^^ 

k 

(where the sum bears on all dyadic cubes of width 2-^ ). For applications, it is important to note that the 
original requirement that the mother wavelet belongs to the Schwartz class can be relaxed to wavelets with 
sufficient regularity and vanishing moments: > |s| if P > 1 and > s > d{2/p — 1) — otherwise [41]. 

When p > 1, Besov spaces are closely related with LP spaces via the following embeddings: 

Ve > 0, ^ 

In particular, the condition rjxip) > 0, used as a practical sharp condition implying that data locally belong 
to LP, is equivalent to requiring that X locally belongs to Bp’°° for an £ > 0. The requirement X G Lp can 
thus systematically be replaced by A G Bp'°° in all results reported so far, with no change in any of them, 
yet with the benefit that it readily extends down to p > 0. The definition of p-leaders given by (21) remains 
unchanged, and p-exponents for 0 < p < 1 can now directly be dehned through (27). 

Admissible distribution. From now on, we therefore assume that there exists po > 0 such that X 
locally belongs to for an e > 0 (i.e., r]x{po) > 0), so that p-exponents can be dehned, at least, for 

0 < p < pq. When this condition is fulhlled, X is said to be an admissible tempered distribution (see Section 
3.5.2 for an example of a random distribution which is not admissible). The extension of Theorem 1, to the 
full range p G (0, -foo], will be the subject of a forthcoming paper. 

Chirps with negative regularity. To illustrate the relevance of using p-exponents with 0 < p < 1, let 
us consider the example of chirp functions (cf. Section 2.4.2). When a < —1, Ca,i 3 does not belong to Lj^^. 
However, Ca^p locally belongs to Bp’°° when p < —l/a and hp{xo) = a. 

3.5. Illustrations, examples and counter-examples 

3.5.1. Illustrations and examples 

The power-law behaviors underlying the wavelet characterization of p-exponents (cf. (27)), and thus their 
p-leader based estimation, are illustrated on the reference examples dehned in Section 2.4, and compared 
with an estimation relying on a direct use of the dehnition of the TP regularity (cf. Section 2.2 and (5) and 
(6)). Fig. 2 shows cusp singularities (cf. (8)) with positive and negative exponents. Fig. 3 illustrates the 
analysis of chirps (cf. (10)), with positive and negative exponents. Fig. 4 addresses the analysis of lacunary 
combs Fa^~/ (cf. (12)), with positive and negative exponents. 

3.5.2. A counter example: white Gaussian noise 

Despite its being very broad, the class of admissible distributions does not contain all tempered distribu¬ 
tions. Notably, the practically very natural example of univariate independent identically distributed (HD) 
Gaussian random variables, also referred to as white Gaussian noise (wGn), satishes = —1/2. However, 
its wavelet coefficients read 

Cj,k = 2 Xj,kj 

where the Xj,k are HD zero-mean Gaussian random variables. This implies that the wavelet scaling function 
reads r]x{p) = —p/2 < 0, for all p > 0, and therefore that p-exponents are not htted to analyze WGN. 
The same holds for fractional Gaussian noise (fGn), the increment process of fractional Brownian motion, 
the only Gaussian self-similar process. Whatever the value of the self-similarity parameter 0 < H < 1, 
Vxip) = —pH < 0 for all p > 0. 
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Figure 2: p-exponents of cusps. Left column: Cusps with positive and negative p-exponent. Second 
column: Estimates of hp{xo) using definition (5-6) (blue circles) and p-leaders (red discs), theoretical 
exponent (black solid) and limit (grey). Third and fourth columns: log 2 Ti^^(a, xq) and log 2 £^^^( 0 , xq) 
as a function of log 2 (a). 



Figure 3: p-exponents of chirps. Left column: Chirps with positive and negative p-exponent. Second 
column: Estimates of hp{xo) using t!^'^ dehnition (5-6) (blue circles) and p-leaders (red discs), theoretical 
exponent (black solid) and limit (grey). Third and fourth columns: log 2 TQ^^(a, xq) and log 2 xq) 

as a function of log 2 a. 
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Figure 4: Lacunary comb Fa,-y{x). Two functions Fa^j{x) with positive (top row) and negative (bottom 
row) p-exponents (left column). Second column: Estimates of hp{xo) using t!^'^ definition (5-6) (blue circles) 
and p-leaders (red discs). Third and fourth columns: log 2 (a, xq) and log 2 (a, Xq) as a function of 
log 2 a. The black solid line indicates the theoretical p-exponent, the grey solid line the limit for the singularity 
to be in 


4. Pointwise singularity classification 

4.I. Motivation 

Traditionally, singular behaviors are categorized using two classes, non oscillating versus oscillating sin¬ 
gularities, cusps and chirps (as defined in Section 2.4) constituting the reference models for each class 
respectively. Deciding whether real-world data contain oscillating singularities or not is of practical impor¬ 
tance, as it may change the analysis of the underlying physical and biological mechanisms at work behind 
data [42, 43]. This is notably the case in hydrodynamic turbulence modeling, where the presence/absence 
of oscillating singularities may permit to validate/falsify various vortex stretching mechanisms [44]. Such a 
classification currently relies on the use of two exponents, the Holder and the oscillation exponents, and is 
best illustrated using the concept of integration (cf. e.g., [26]). 

Let us compare the ID chirp function (cf. (10)) to its primitive For the former, the p- 

exponents at Xq read hp{xo) = a (for p € (0,-Foo] if a > 0 and p G (0,—1/a) if a < 0) and are hence 
independent of /? and p. For the latter, a simple integration by parts yields 

V^^p{x) = -cos + C + 0{\x- xo|“+2^+"), 

thus showing that the increase of any p-exponent after integration reads j3 -\- 1. This is in contradistinction 
with the cusp case (8), where the increase after integration is exactly 1. This simple computation indicates 
that cusps and chirps can be discriminated by integration and paves the way towards the definition of an 
oscillation exponent (3 (cf. e.g., [26] for the case of Holder exponents). 

However, both cusps and chirps are characterized by p-exponents that do not depend on p and can hence 
not be considered different with respect to p-exponents, but instead depart from the lacunary comb example 
(cf. Section 2.4.3). A general interpretation of these two different behaviors is provided in [43]. 

This naturally leads to define the new class of p-invariant singularities. 

Definition 2. Let X : M. be an admissible distribution. The singularity at xq is said to be p-invariant 

if and only if 

Vp G (0,po(xo)), hp{xo) = h. 
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In other words, the function p —>■ hp(xo) is constant in the interval (O,po(a;o))- 

The remainder of the section aims to consider the extent to which the flexible setting of p-exponents 
permits to enrich the classification of singular behaviors, as illustrated in Fig. 5. 


4-2. Canonical vs. oscillating singularities 
4-2.1. Fractional integration 

While examples of the previous section were ID functions, the multidimensional setting, a; G with 
d > 2, is considered, thus requiring to replace the concept of primitive function with that of fractional 
integral: 

Definition 3. Let X be a tempered distribution defined on The fractional integral of order s of X, 
denoted by X^~^\ is defined as the convolution operator (/d — Equivalently, in the Fourier domain, 
it corresponds to the multiplication by the function (1 + The p-exponent of X^~‘^'> at xq is well 

defined on condition that rjxijf) > —sp and is denoted by hp^sixo). 

The condition rjxip) > —sp, sufficient to insure that X‘-~^'> locally belongs to when p > 1 (or to Bp’°° for 
an e > 0 when p < 1), follows from the Besov space interpretation of the wavelet scaling function, see [35]. 
Generalizing (27), the exponent hp^sixo) can be characterized by 


with (p, s)-leaders defined as: 


hp,s{xo) = hmmf -- ■ ■ ■ 

J^-oo log(2-?) 


2“-l 


1/p 


= 4^’"^ = I Y (2'^" l^v 1)^ 

KT<3, A'C3A i=l 


(29) 


(30) 


Essentially, this means that, as regards pointwise regularity, taking a fractional integral of order s is equiv¬ 
alent to multiplying the wavelet coefficients by 2-1®, see [31, 40] for a mathematical justification of this 
heuristic. 


4-2.2. Canonical singularity 

The notion of canonical singularity is defined as follows: 

Definition 4. Let p and hp be such that p < Po{xq) and hp > —d/p. An admissible distribution X : —>■ K 

has a canonical singularity of exponent hp at Xg if Vxip) > 0 and if there exists s > 0, such that 

hp^s{xo) = hp + s. (31) 

A similar notion had previously been considered only in the Holder case, i.e., when p = -|-oo, see [45]. The 
following result shows that canonical singularities are p-invariant singularities (thus motivating the choice 
of the terminology). 

Theorem 2. Let X : K®* — > M 6e an admissible distribution with a canonical singularity at xq. Then 

VpG(0,po], Vs > 0, hp^s{xo) = hp + s. (32) 

In particular, p hp{xg) is constant in the interval (0,po(xo)). 

The proof of Theorem 2 is detailed in Appendix C. Note that it does not rely on the assumption that 
Po(xo) = +00, and its conclusion also applies for 0 < p < 1. 

Example 4: Self-similar distributions. Besides cusps, further examples of canonical singularities 
are supplied by deterministic self-similar functions or distributions, which are defined as follows. We only 
consider the one-variable case here. 
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Definition 5. Let a > 1 and a S M. ^ nonvanishing tempered distribution X defined on M is self-similar 
of scaling ratio a and exponent a at xq if the two distributions a°'X{x — xq) and X{a{x — xg)) coincide. 


The following result shows that self-similar functions supply simple examples of canonical singularities. 


Proposition 1. Let p > 1 and X € be a self-similar function at xg of scaling ratio a > 1 and 

exponent a. If X differs from a polynomial in the neighbourhood of xg, then p > —1/a, X has a canonical 
singularity at xg of exponent a, and it can be written under the form 


f lx - xo|“a;+(log(a; - xo)) if x>0 

\ jx - xoi“a;_(log(-(a; - xo))) if x < 0, 


(33) 


where a;_|_ and uj- are log(a) periodic functions in L^. 


The proof of Proposition 1 is given in Appendix D. 

Simple examples of self-similar functions or distributions are provided by the cusps Ca when a ^ Z 
(defined by (8) for a > —1 and in [43, Section 3.4] for a < —1), in which case any a > 1 will fit. Note that, 
if a > — 1, these fall in the case covered by Proposition 1 with uj+ and oj- set to constants. Several examples 
of cusps are plotted and analyzed in Fig. 2. Another important example is provided by singularities with 
complex exponent [46] 

|x-xor+*^ = |x-xorexp(-j/31og|x-xo|). 

Other examples of self-similar distributions which are not functions are supplied by the measures 

for a > 1 and a > 0. 


Note that, if A is a self-similar distribution of scaling ratio a and exponent a, then its derivative X' (in the 
sense of distributions) is self-similar of scaling ratio a and exponent a — 1; thererefore, by derivation, we can 
deduce from the previous examples self-similar distributions of arbitrary negative order. 


4-3. Oscillating singularities: Lacunary versus balanced 
4-3.1. Oscillating singularities 

Definition 6. Singularities that are not canonical are called oscillating singularities. 

The classical examples of oscillating singularities are provided by the chirps Ca^p in (10). Their p- 
exponents are given by hp(xg) = a (for any p G (0,-|-oo] if a > 0, and for p G (0,—1/a) if a < 0). The 
functions Ca,p therefore supply examples of singularities which are not canonical (they are oscillating) yet 
are p-invariant. 


4-3.2. Lacunary versus balanced singularities 

However, in general, oscillating singularities need not be p-invariant and this large class of singularities 
can be further refined by the use of p-exponents. Let us consider again the function defined in Section 
2.2 and (12). It has an oscillating singularity at xq: Indeed, the primitive of Fa^y has exponent a -I- 7 at 
xo = 0. This directly follows from the fact that its integral on the interval [2“*, 2“* -|- 2“'*'^] is and 

that Fa^y vanishes elsewhere. The major virtues of the examples as well as of the more general case 
developed in Appendix B, consists in showing that observing p-exponents that differ is deeply tied to the 
notion of spatial/temporal lacunarity: Indeed, the support of the function Fq, y, is a set which is “scarce” at 
Xo (cf. [43] for the technical formulation of scarcity). 

This leads to refine the classification of oscillating singularities. 

Definition 7. Let X be an admissible distribution, with an oscillating singularity at Xg. This oscillating 
singularity is balanced if the function p 1— hp(xg) is constant in its interval of definition, and is lacunary 
otherwise. 


The p-exponent based classification of singularities provided by Definitions 2, 4 6 and 7 is schematically 
illustrated in Fig. 5. 
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Figure 5: Classification of singularities. 


J^.^. Negative regularity 

4.4- 1. Motivation 

Let us now return to the original motivation of analyzing negative regularity singularities. As mentioned 
in the introduction, the classical strategy to handle negative singularities is to compute h™'" (using (23)) 
and, whenever /i™™ < 0 to perform a fractional integration of order s > —Indeed, the function 
satisfies + s, which is positive, so that the Holder exponent can be used in order to measure 

its pointwise regularity. A major advantage of this approach stems from its not relying from any a priori 
assumption. A significant limitation stems from the fact that (fractional) integration may mask some of the 
regularities existing in original data. Let us examine how the alternative p-exponent approach may overcome 
such limitations. 

4.4- 2. Example 5: Cusp plus Chirp 

We now construct a simple toy-example that will show why pointwise regularity information can be lost 
in (fractional) integration, when data contain oscillating singularities and negative regularity, cf. [47]. Let 

Xi{x) = \x- xoP, X 2 {x) = \x- xo\°‘ sin ( .-^ ) , and X{x) = Xi{x) + X 2 {x), 

\\x-Xo\P J 

with a < 7 < < 0. The condition a < 7 implies that the dominating singularity at Xq is supplied by A 2 . 

The exponents of Xi and A 2 are respectively 7 and 7 ^, thus yielding /i™'" = 7 for X. Interestingly, 
^mm jg thus larger than a and hence no longer yields a lower bound of possible pointwise exponents, as is 
the case when only singularities with positive regularity are present. A Holder regularity based analysis of 
X requires performing a fractional integration of order s > — 7 . Then, the Holder exponent of X\ at Xg is 
shifted to 7 -F s, and to a -F s(l -F /3) for X 2 , which is strictly larger than 7 -F s. Therefore, the oscillating 
part X 2 ' is dominated by the cusp part A] , and the information associated to the initially dominant 
oscillating singularity cannot be recovered. Note that in this example, p-exponents are equally affected by 
this limitation caused by the fractional integration. 

In contrast, we can use p-exponents with p < —1/a without a prior fractional integration. The dominant 
singularity (chirp) is a and p-exponents enable its correct estimation. This example is illustrated numerically 
in Fig. 6 . To conclude, in the presence of oscillating singularities, using p-exponents should be preferred 
to the recourse to fractional integration, which induces smoothing and may make some of the singularities 
actually existing in data invisible. 

5. Conclusions and Perpectives 

As an alternative to the traditional Holder exponent, the present contribution has put forward the use 
of a collection of new indices for pointwise, or local, regularity quantification: p-exponents, p > 0, allowing 
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sum of cusp and chirp singularities 


fractional integral of order s = 0.3 


Chirp+Cusp 


a=-0.3 

(3=0.4 

Y=-0.25 






Figure 6: Chirp and cusp. Sum of cusp and chirp singularities (top row). The p-leader based estimation 
of p-exponent (bottom row) for original functions (left column) and functions after a fractional integration of 
order s = 0.3: theoretical value for hp{xo) for cusp (blue dashed line) and chirp (black solid line), estimates 
(red disks) and limit (solid grey). The value for p = oo corresponds to the Holder exponent. Only small 
values of p allow to recover correctly the order of the dominant singularity: a = —0.3 


notably the practical use of the notion of negative regularity, of significant importance for real-world data 
and applications. A corresponding wavelet framework enabling the practical estimation of p-exponents, the 
p-leaders, has been theoretically studied and numerically illustrated at work. 

The construction and analysis of p-exponents and p-leaders has yielded additional comprehensive under¬ 
standing on local regularity analysis. First, it has shown that the choice of a particular index to quantify 
local regularities forces the selection of the relevant multiscale quantities for theoretical and practical anal¬ 
yses (while it is often misleadingly understood that the choices of the regularity index and of the multiscale 
quantities — increments, wavelet coefficients — can be independently achieved): While the Holder exponent 
is intimately associated to oscillations and wavelet leaders, p-exponents require the use of p-leaders. Second, 
the present work has shown that all regularity indices are not equivalent in general; p-exponents hence 
permit to enrich the classification of singularities in terms of p-invariant versus non p-invariant, canonical 
versus oscillating, and oscillating balanced versus oscillating lacunary behaviors. 

The present contribution focused on the definitions and properties of p-exponents and p-leaders. It has 
hence been illustrated on isolated singularities for the sake of pedagogical exposition. The general aim of 
the present work however is multifractal analysis, aiming to study collections of intertwined singularities. 
The construction of a p-leader based multifractal formalism is the topic of a companion article [38] , which 
illustrates the additional benefits of p-leaders even in situation where p-exponents and Holder exponent 
coincide. This new formalism will also permit to show how p-leaders connect, generalize and outperform 
theoretically and practically Multifractal Detrended Fluctuation Analysis^ an earlier multifractal formalism 
of large popularity notably in biomedical applications (cf. [3] for the seminal definition). 

Matlab routines designed by ourselves implementing p-leaders and the estimation procedures for p- 
exponents will be made publicly available to the research community at the time of publication of the 
present article. 
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Appendix A. Proof of the conditions satisfied by the p-exponent function (first part of Theo¬ 
rem 1) 


Proof: If A G then 


/ B{0,a 


\X{u + xo)\^du < C 


so that the left hand-side of (5) is bounded by thus X G T^^i/pixo) and the first point of Theorem 1 

holds. 

Assume that X G TP{xq) for a p > 1 and let q be such that 1 < q < p. We denote by P the Taylor 
polynomial of X for the TP{xo) regularity condition. Let B = B{xo,r), p' = p/q and q' the conjugate 
exponent of p', i.e. q' satisfies ^ ^ = 1- Using Holder’s inequality, 


f \X{x) — P{x — xo)\‘^ dx = f \X{x) — P{x — xo)\‘‘ 1b{x) dx 
J B J B 




1/p' 


J \X{x) - P{x - xo)\'^P dxj (Ib(x))'^ dx 

J \X{x) — P{x — xo)\^ dx^ (C'r'^) 


1/9' 


>1/9' 


_ QpCiq+d 


It follows that X G T^{xo) and, up to the order hp{xo), the same Taylor polynomial can be used for T^{xo) 
as for TP{xo). Hence the second point of Theorem 1 holds, and also the assertion concerning the Taylor 
polynomial in Section 2.2. 

In order to prove the concavity of the function s —>■ hi/s{xo), we assume that X G TP{xo) n Tp{xo). We 
pick 7 G (0,1) and define s by 

I 7 1 — 7 

^ = —I-; 

s p q 

we have to prove that X G Tg{xo) where 


^ = 7a -I- (1 — 7)/?. 

We pick for P the Taylor polynomial of largest possible degree of X. Denote by || X ||p the norm of X 
on B{xo,r); the assumptions imply that 

\\X-P ||p< Cr^^+d/p and || A - P ||g< 

By interpolation, it follows that 

II X — P ||s< (7pT'(“+'=*/P) + (l“')')(/5+‘^/9) = (7p7a + (l-7)/5+d/s 

SO that A G Tg{xo). Hence the concavity of s — ?► ^i/s(a^o) follows. 
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Appendix B. Proof of the optimality of Theorem 1: Construction of functions with arbitrary 
p-exponents 

We now show the optimality of Theorem 1 by showing that the function p —>■ hp{xo) can be any function 
satisfying the conclusions of Theorem 1. The proof consists of a constructive example. 

For the sake of simplicity, we do the construction in dimension d = 1; the reader will easily check that it 
extends to arbitrary dimensions. 

The functions that we will consider are parametrized by two sequences 9{l) and uj(l) which satisfy the 
following conditions: 

1. ^ ^ “t-oo 

OO 

2. < OO 
1=1 

The function is defined as follows (we take here a;o = 0) : 

if a: e [2-', 2-' + 2-‘^(')] for an integer I > 0, then Fe^^{x) = 

otherwise Fg^^{x) = 0. 


Note that the first condition that we imposed implies a lacunarity on the construction: As a: —^ 0, Fg^i^ 
vanishes on a larger and larger proportion of points. The second condition implies that Fg.tj belongs to 
in the neighbourhood of 0. 

First, it is clear that Fg^^^ is bounded in the neighbourhood of Xg if and only if 

3(7 € K such that VI > 0, 0(1) > C. (B.l) 

Lemma 1. If (B.l) holds, then the Holder exponent of Fg^^ at 0 is given by 

hg,uj{0) = liminf (B.2) 

I —^-|-00 i 

Proof: Indeed, we first note that Fg is continuous at 0. Remark also that for a > 0 the local maxima 
of the function Fg^,^(x)/x'^ are obtained at the points 2~K Thus if we pick P = 0 in (3), then computing 
the Holder exponent amounts to find the supremum of the a such that one can find (7 > 0 such that for all 
I > 0 < C. This is easy to check that this supremum is exactly given by (B.2). 

We still have to check that the estimate obtained by taking P = 0 is the best possible. As regards the 
constant term, we separate two cases: if 0 (1) does not tend to +oo then it follows that /ie,tj(0) = 0; and if 
0(1) —>■ +00, the constant term necessarily vanishes. For higher order terms, one argues by induction on the 
valuation of P, noticing that Pg^,^ vanishes on the intervals [2“^ + 2“*+^], which implies that, on such 

interval, the choice of non-vanishing terms for P would lead to worse estimates. 


We now come back to the proof of the second part of Theorem 1 and we no longer make the assumption 
that (B.l) holds. In order to estimate the p-exponent of Pg^^^ at 0, we start by considering the quantity 


1 

a 


\Fg^^(x)\Pdx 


(B.3) 


(we thus take P 
then 


0 in the definition of the Tf) regularity). First assume that 2 *“ + 2 < a <2 


1 

a 


\Fe^^(x)f‘dx 



(B.4) 
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Therefore, the critical value p(0) is given by 

p(0) = sup{p : (B.4) converges}. 

We will now make explicit choices for the sequences 9{l) and uj{l). We pick the variable s = 1/p, so that 
the conditions satisfied by the function 

p{,) = 4 , Q) 

are that it is a concave increasing function defined on (l/p(0), 1] and satisfying 


P{s) > 


1 

p(0)' 


Concavity implies that p can be obtained as the infimum of a countable family of affine functions 

Pn{s) = a„s + bn 


which all satisfy 


Vs,p„(s) > p{s). 


(B.5) 


We can also assume that there exists a dense sequence Sn such that Pn{sn) = p{sn)- We now pick the 
functions 9 and uj as follows: Note that any integer I > 1 can be written in a unique way under the form 
I = 2"(2fc + 1), with n > 0 and fc > 0; then 


f ^(0 “ + 1 )^ 

1 9(1)= bnl. 


The result follows from the 


fact that partial sums of the righthand side of (B.4) satisfy 


1 1 

_ _ o-‘^9)o-p^9) < 

l = lo 


1 

a 



— (a„ + l+p6„)/ 


which, using (B.5) and p > 1, is bounded by 


(B.6) 


m - 

C2’‘° _2~^PP^'^/p)+'^)^ 
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l = lo 


Using that the exponent is strictly positive (since 4> is defined on [l,p(0)]), this sum is bounded by 


(^2'o2-(pp(i/p)+i)^o < (j ( qP(i/p) 


p 


(B.7) 


which is the required upper bound. 

If we now take a in the interval [2“*“, integral is bounded by the value at 2~^° 

which has already been estimated, so that (14) still holds. We still have to check that the estimate obtained 
by taking P = 0 is best possible. As in the Holder case, one argues again by induction on the valuation 
of P, noticing that Fg ^ vanishes on the interval [2“^° + 2 ~‘^^''°\2~''°], so that the order of magnitude of 
I\Fe, c^{x) — P{x)\Pdx on such intervals will be given by the integral, on this interval, of the first non-vanishing 
term of P. This remark immediately yields that, indeed, the choice P = 0 is best possible. 


The lower bound is obtained by noticing that the first term of the sum in (B.6) is larger than 

1 1 




2“(“ii + l+Pfenl^O 
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which for p = pn ='■ ^/Sn and a = 2 + 2 is larger that 


1 ^_O-PnP(l/Pr,)i0. 

the lower bound is thus sharp at the points p„. Since they form a dense set, and since p is continuous, (7) 
is proved. 

We have therefore obtained that the p-exponent of Fo^^i at the origin is 

/ip(O) = p{l/p), (B.8) 


and Theorem 1 is proved. 


Appendix C. Proof of Theorem 2 


Proof: We start by introducing a useful notation. Let p S (0, +oo); we will say that the sequence 


y Aj (ico) J 


i<o 


satisfies 




j<0 


2^3 


at xq 


if the two following conditions hold: 

Ve > 0, for j small enough, ^ 


Ve > 0, 3jn —> —00 
Note that these two conditions can be rewritten as 


f(p) 

Aj„(a:o) — 


> 2 (h+e)j„ 


lim inf 

j^-oo 


■°8 «*(..)) 

log(2J) 


= h: 


(C.l) 

(C.2) 

(C.3) 


and, since r]x(p) > 0, the wavelet characterization of p-exponents given by (27) exactly means that (C.3) is 
equivalent to hp(xo) = h. It follows from (29) that X has a canonical singularity of exponent hpg at xq if 
and only if 

47io) "" at xo and 3s > 0: 47(£) ~ at Xq. 

Theorem 2 will be a direct consequence of the following result. 

Proposition 2. An admissible distribution X has a canonical singularity of exponent hpg at Xg if and only 
if (C.l) holds for p = po and if, for any e > 0, there exists a sequence of dyadic cubes A„ of scales jn —>■ —oo 
such that 

dist{xo, Xn) < 

(C.4) 

dA„ > 2^^P0+C3r, 

where d\ is the quantity defined in (22). 

Proof of Proposition 2: Let us assume that X has a canonical singularity of exponent hpg at Xq. We 
will first prove that 


l/po 


Ve > 0, 3jr. 


—oo 


E 

\\' C 3 \ j { xo ), jn>j'>(l+e)3n 




1 - 3 ') 




(C.5) 
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We prove it by contradiction. Indeed, if it were no true, then there would exist e > 0 such that. 


l/po 


Vj<0 ( 

yA'C3Aj(a;o). j>j'>il+e)j 

Let s > 0; consider the corresponding (po, s)-leader (30): 

V A'C3A„ (xo) 


< 2^^ PQ +e)j ^ 


1/Po 


(C.6) 


(C.7) 


we split the sum into two parts, depending whether j' > (1 + e)j or j' < (1 + e)j. 

It follows from (C.6) that, in (C.7), the term corresponding to j' > (1 + e)j is bounded by 2(^po+^+®)J. 
Consider now the term corresponding to j' < (1 + e)j; since 2®-^ < it follows that it is bounded by 

2(^Po+®(i+E))i_ Therefore 

hpo,sixo) > hpg + s + e(min(l, s)), 

which contradicts the fact that xq is a canonical singularity. Therefore (C.5) holds. 


We now prove that (C.5) implies that (C.4) holds. Let e' > 0. We first estimate the number of coefficients 
on which the sum in (C.5) bears using e'. At the first generation j' = j„ — 1, there are 2 ‘^ subcubes of (xq); 
at each generation, each subcube of the previous one is split into 2 ‘^ cubes; we go down to the generation 
j' = [(1 + e')jn] = jn + [s'3n]^ where there are (2‘^)[“®'^"1 cubes. Therefore, the total number of cubes 
considered is bounded by 3 • (2^^)!“® 

It follows that one of the terms of the sum in (C.5) satisfies 

which, using that j„ > j' > jn + s 'implies that 

dy > C2 ^^o3'2^'^'^'. 

Furthermore, the conditions X' C 3Aj^(a;o) and jn > j' > jn + s' jn together imply that the cube X' satisfies 
the first condition of (C.4). 


We now prove the converse part in Proposition 2. Assume that (C.l) and (C.4) hold. We denote by A 
the smallest cube of the form Aj(xo) such that 3A contains A'. It follows from (C.4) that the scale j of A 
satisfies j'{I — s') < j < j'. Therefore the corresponding pg-leader satisfies 

> 2('»P0+®')/2'^®'^ > (C.8) 

and, similarly, the corresponding (po, s)-leader satisfies 

Remark that any e > 0 can be written Ce', thus hp^^^ixo) < hp^ + s. Note that (C.l) implies that 
hpoixo) > hp„, and the lower bound hp^^g > hp^ + s then follows from general results on fractional integra¬ 
tion, see [31, 40]. We have thus obtained that X has a canonical singularity of exponent hp^ at a:o- 


J-/P0 




(C.9) 


We now turn to the proof of Theorem 2. First, we note that the derivation of (C.8) also holds for p < po, 
and it follows that 

Vp < Po, hp{xo)<hpg. 

Since p i—>■ hp(xo) is decreasing, it follows that Vp < po, hp{xo) = hp^. Since (C.8) also holds for p < po, by 
the same arguments as above, we deduce that 

Vp < Po, Vs > 0, 
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hp^s[x(}) — ^pq 4 “ 


Appendix D. Proof of Proposition 1 


Proof: We can assume that xq = 0. We first prove the second part of the proposition; let us define the 
functions uj+ and tO- by (33); this implies that they belong to on [0,loga]. The self-similarity condition 
exactly expresses the fact that they are logo periodic functions. The condition a > —1/p then follows from 
the condition that X G in the neighbourhood of 0. 

We now compute the p-exponent at Xq. We first assume that the Taylor polynomial of X vanishes at 0. 
For r = 2“*^, 

\X\Pdx = 



\x\^P{uj+{\ogx)Ydx 



\xY^{uj+{\ogx)Ydx 


oo „2 

= 2 ^ ^ 2 -^ / 2-‘^P^\uYP{uj+{\ogu)Ydu 

0=J+1 “'l 

which is finite if and only if op > —1, in which case its value is C ■ The quantities are of the same 

order of magnitude if r lies between two powers of type 2“'^, and the computation is the same for a; < 0. 
Therefore the p-exponent of X at 0 is a. 

We now check that the Taylor polynomial of X vanishes at 0. First, if a is not an integer, this is clear, 
since a polynomial of degree less than a would bring a higher order contribution in the above computation. 
For the same reason, if a is an integer, then the polynomial necessarily is a monomial of order a. Now, if 
P{x) = then we note that X — P also is self-similar of order a, so that the above computation still 

applies, and the integral can be of a smaller order of magnitude only if X coincides with P on both sides of 
the origin, which we excluded in the assumptions. 

Let us now prove that X has a canonical singularity at the origin. We note that, through an integration 
by parts, one easily checks that, if X is a self-similar function of scaling ratio a and exponent a at xq, then 
its primitive which vanishes at 0 is self-similar of scaling ratio a and exponent a -I- 1. It follows from the 
previous computations that its p-exponent is a -t- 1, so that (31) holds with s = 1, hence Proposition 1 is 
proved. 
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